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ABSTRACT 

We construct infrared star count models of the Galaxy applicable at faint 
magnitudes based on the models and luminosity functions obtained at V. There are 



^ ■ two types of models derived in our paper. First we derive infrared versions of the 

disk and halo models by Gould et al. who obtained the local luminosity functions 
i and functional forms of the models based on HST star counts (HST based models). 

^ I Second we construct a double exponential disk model based on the luminosity function 

' of the nearby stars by Reid and Gizis and a halo model based on the nearby subdwarf 

' luminosity function of Dahn et al. (Nearby star LF based models). In addition to 

00 ' dwarfs and subdwarfs included in original models, we also take into account L dwarfs, 

■ T dwarfs, and white dwarfs of both disk and halo. As a test case of the models, we 
Q . analyze the infrared imaging data at J and K' obtained during the Subaru Deep Field 
^ ! (SDF) survey to study stellar objects. Out of some 350 objects, 14 stellar objects are 

Qh| selected in a 2' x 2' field based on a morphological criterion applied to the J band 

O I image. Both completeness and contamination associated with the selection criterion 

' are evaluated by simulations. The J band image is 57% complete at J = 24 and the 

■ number of contaminants is estimated to be negligible. The prediction of the HST based 
^ ■ models agrees with the observed count at — O.80" and that of the nearby star LF based 

^ . models also agrees with the observations at +1.0a. The observed count is between 

the predictions of the two types of models that have contradictory local luminosity 
functions. With our limited statistics, the observing data do not favor a particular 
type of models. Infrared star count models we have obtained predict that the Next 
Generation Space Telescope will primarily see T dwarfs, M subdwarfs, and old halo 
white dwarfs at faint magnitudes. 
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1. Introduction 

Deep infrared images have been obtained so far for the purpose of deep galaxy counts (e.g. 



Bershady et al. (1998) ). It is indeed the case that the majority of the objects detected in those 
images are galaxies and stellar objects are the minority. Infrared images of the Subaru Deep 
Field (SDF), which we are going to analyze in this paper, were obtained also for the purpose of 



a galaxy count analysis ( Maihara et al. (2000)[ ). However as we show later, these images contain 



interesting information about faint stellar /substellar populations, since they reach the magnitude 
limits which no stellar astronomer has ever experienced. 

The SDF images were obtained at J and K' . Once stellar objects are selected from a catalog of 
detected sources by a criterion on the morphology of individual sources, only remaining information 
is color-magnitude information. As an aid to interpret the color-magnitude information, we 
extensively utilize Galaxy model predictions on the stellar/substellar populations. Infrared star 
count models at faint magnitudes have not been constructed previously. We generate models for 
J based on those constructed for V and color-magnitude relations for visual and infrared. We 
also take into account faint stellar/substellar populations such as L dwarfs, T dwarfs, and white 
dwarfs of the Galactic disk and halo. We compare the model predictions with observations and 
examine the consistency. We also utilize these models to predict the Galaxy observed by the Next 
Generation Space Telescope (NGST). 

In the past, infrared images as deep as -ftT = 19 were analyzed from the point of view of 
stellar astronomy. Hu et al. (1994)| studied the deep imaging data obtained for a deep galaxy 



count and concluded that very little of dark matter in the Galactic halo can be made of low-mass 
hydrogen-burning stars. In our small field of view, there are only 14 stellar objects. However 
analyses of a small number of stellar objects have been already done in the Hubble Deep Field 
( [Elson et al. (1996) , [Mendez et al. (1996) ). Using different criteria, Elson et al. found 59 



unresolved objects while Mendez et al. selected only 14 stellar objects. Elson et al. counted 
only nine objects as candidates for low mass halo stars and ruled out a steeply rising luminosity 
function fainter than My = 12. Mendez et al. found that the luminosity function for halo objects 
must be a factor of 2 smaller than that predicted by an extrapolation of the solar neighborhood 
luminosity function for disk stars. 

Star count models are constructed in §2. The reduction and analysis of the SDF data are 
described in §3. Model predictions are compared with observations in §4. Individual objects in 
SDF are discussed in §5 and the Galaxy seen by NGST is predicted in §6. 



2. Star Count Models 



Infrared star count models applicable to very faint magnitudes do not exist in the literature. 
So we construct them, starting from the models obtained in the visible, using color-magnitude 
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relations. We also use a K band luminosity function for the disk as a starting point when 
available. We obtain color-magnitude, magnitude-magnitude relations empirically from published 
photometric data when necessary. We first construct models by Gould et al. who used HST 
star count data. Their local luminosity functions (LLFs) of the disk and halo obtained by fitting 
distant stars are inconsistent with the LLFs of nearby stars. We also construct models different 
from those of Gould et al., using the nearby star LLFs. 



2.1. Model A: Disk model by Gould et al. 

We first construct an infrared version of the disk model by (Gould et al. (1997)) who 



fitted the disk structure parameters and the LLF to star count data from HST observations of 
high-latitude fields. The distribution of stars is modeled as a function of Galactic position and 
absolute magnitude by 

^{Mv,z,R) = ^Mv)v{z)eM-^^-^). (1) 

where ^{My) is the LLF, {R,z) is the Galactic position in cylindrical coordinates, Rq = 
8 kpc is the Galactocentric distance, and H = 2.92 kpc is the disk scale length. The vertical 
distribution function is given by 

z^(z) = (l-/3)sech2-^+/3exp^, (2) 
hi h2 

where hi = 320 pc, /i2 = 643 pc, and /? = 21.6%. 

The local disk/thick disk ratio, /3, adopted by Gould et al. is more than a factor of two 
higher than any other estimates and this discrepancy is reflected in the very small scale height, 
/i2- This small scale height results in a smaller star count estimate at faint magnitudes. The LF 
is given for main-sequence stars with 8.0 < My < 18. Since we are interested in star counts at 
faint magnitudes, this range of the LF is appropriate. The infrared color-magnitude relation is 
not available for thick disk stars and we use the color-magnitude information of thin disk stars in 
Leggett (1992) for the entire disk stars to calculate star counts in the infrared. Reid et al. (1997) 



have found that the metallicity of the stars do not change significantly out to the scale height of 
2 kpc. This may be a support for the use of the thin-disk color-magnitude relation for the entire 
disk stars. 
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2.2. Model B: Halo model by Gould et al. 



We also construct an infrared version of the halo model by Gould et al. (1998) who fitted the 
halo structure parameters and the LLF to HST star count data. The distribution of halo stars is 
given by 



i?2 + {z/cf 



Rl 



-1/2 



(3) 



where ^hi^v) is the LLF, c = 0.82 is the flattening parameter, and / = 3.13 is the power. 
There is no infrared color-magnitude relation available for halo stars. To calculate star counts 
at /, J, K and J — K color, following empirical relations were obtained by fitting extreme M 
subdwarf data in Leggett et al. (1998)| ; 



Mi = 0.745Mi/ + 1.122, 



(4) 



Mj = 0.638Mi/ + 1.249, 



(5) 



Mk = 0.649Mi/ + 0.503, 



(6) 



J -K 



-O.OllMy + 0.747. 



(7) 



The sample in Leggett et al. also includes a small number of M subdwarfs (not extreme M 
subdwarfs), but they are not used in the fitting. 



2.3. Model C: Exponential disk model based on nearby star LF 



We construct a double exponential disk model using the nearby star K band LLF by Reid 
and Gizis (1997)| . The distribution of disk stars is given by 



^{Mk,R, z) = $(Mx) exp 



^)[(l-/3)exp(-M) + ;9exp(-M] 



(8) 



where ^{Mk) is the K band LLF, i?o = 8 kpc is the Galactocentric distance, H = 3.5 kpc is 
the Galactic scale length, hi = 325 pc is the scale height of thin disk stars. For the scale height 
/i2 and the normalization (3 of the thick disk, we consider the combinations of (/i2i/3) = (910 pc, 
(Model CI: puser et al. (1999)| ) and (/i2,/?) = (1500 pc,2%) (Model C2: |Reid and Majewski| 
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(1993)). Using the color-magnitude relation of thin disk stars by Leggett (1992) , the model is 
constructed so that it predicts star counts at V, /, J and K. We note that our model has a 
limitation in the use of the thin disk color-magnitude relation for the thick disk. However the 
effect of metallicity on the J — K color is relatively small and the J band count is not significantly 
affected by this problem. 



2.4. Model D: Halo model based on nearby star LF 



The model D is based on the LLF of nearby stars by Dahn et al. (1995). The density 
distribution of halo stars is given by 



^{Mv,R,z) 



HMv){^) '^'|exp[-10.093(^)V4 + 10.093]} 

X [l 



X 
-Rq 

0.08669/( 



^ ^1/4 

Rq 



(9) 



where x = {R"^ + (z/cf)^/^ and c = 0.85 ( |Reid and Majewski (1993]| ). $(Afy) of Dahn et 
al. is significantly greater than the LLF of Gould et al. (1998)[ Star counts at /, J, and K are 
estimated from the magnitude-magnitude relations eqs (6), (7), and (8). 



2.5. L dwarfs, T dv^arfs and vi^hite dvi^arfs 



Dwarfs later than M, namely L and T dwarfs ( Kirkpatrick et al. (1999)| ) are also included 
in our models. The local number density of L dwarfs is estimated by the 2MASS observations of 

max estimates. The local density of T dwarfs is 



Kirkpatrick et al. to be 7.2x10 ^ pc ^ using 1/V 



estimated by the SDSS observations of ptrauss et al. (1999)| and p?svetanov et al. (2000)| to be 



1.4x10"^ pc"3. This assumes that the survey volume extend out to 10 pc. As for T dwarfs, we 



note that the 2MASS observations by Burgasser et al. (1999) give a lower number density than 
the SDSS observations. For simplicity, L dwarfs are assumed to be a single population with Mj = 
13.0 and J - = 1.5, and T dwarfs with Mj = 15.4 and J - K = -0.1 like Gl 229B ( |Matthew£ 



et al. (1996)). The ages of L and T dwarfs are expected to be order of 1 Gyr and they belong to 
the young disk population. The scale heights of both L and T dwarfs are assumed to be 200 pc or 
the scale height of F stars. 



We then consider white dwarfs. For hot white dwarfs, we compute the LF by [Liebert et al. 



at V. Since there is no good observational color-magnitude relation for hot white dwarfs 
usable to derive an infrared luminosity function, we are satisfied with the V band model to order 
estimate infrared star counts. The luminosity function of cool white dwarfs {M{,oi > 13) has been 
obtained by Leggett (1998) for Mf,oi based on visual and infrared observations. By fitting her 
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observational data, we have obtained following relations among M^oi, My, Mi, Mj, Mk and 
J-K\ 



My = l.l03Mhoi - 1.311, (10) 

Mi = 0.743Mboi + 3.084, (11) 

Mj = 0M2Mboi + 5.593, (12) 

Mk = OAMMboi + 6.041, (13) 

J-K = 0M7Mboi - 0.448. (14) 

Using these relations, we include cool white dwarfs in the disk models. It turned out that 
disk white dwarfs do not have significant contribution at faint magnitudes. 

We also consider old white dwarfs in the halo ( [Hodgkin et al. (2000)| ). These white dwarfs are 
very faint and blue in the infrared. The absolute magnitudes are assumed to be of WD 0346+246 



( Hodgkin et al. (2000) ) and the number density is assumed to be 7x10 pc ( Ibata et al 



2000)1). These old halo white dwarfs turn out to be important in the star count by NGST at the 



faintest magnitudes as will be discussed in §6. 



3. Reduction and Analysis of the SDF data 

Details of the observations of the Subaru Deep Field (SDF) are described in [Maihara et al.| 
'2000)1 . The SDF is a 2' x 2' field centered on {a,6) = (13''24'"21.3^ 27°29'23") (J2000) or (l,b) 



= (35.°7, 82.°2). For the purpose of star-galaxy separation, exposure frames taken only under 
good seeing conditions were coadded to form final images. The total integration time of the 
coadded frames at J is 4.97 h out of 12 h of observations and that at K' is 4.13 h out of 10 h of 
observations. The FWHMs of the resultant PSFs are 0."39 and 0."29 respectively at J and K'. 



A source catalog was created for each band by the SExtractor ( Bertin and Arnouts (1996)| ). 
For each source, its position, magnitude, FWHM of the image, and 'stellarity' were recorded. 
The stellarity is a parameter defined in the SExtractor which quantifies how stellar the object 
morphology is. It is unity for a completely stellar object and zero for a completely non-stellar 
object. The stellarity is an output parameter of a neural network. The purpose of the use of the 
neural network is to provide an optimal transformation from the parameter space defined by a set 
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of observables describing the object, to the one dimensional space of stellarity. Input parameters 
are eight isophotal areas, the peak intensity, and the seeing FWHM. 

Star-galaxy separation is a critical process in the analysis of a deep image. Its difficulty in a 
deep optical star count has been mentioned by Reid et al. (1996), who analyzed visible images 



obtained on the Keck telescope under 0."5 seeing. They conclude that the decreasing average 
angular size of galaxies with fainter magnitude limits star-galaxy separation to R < 25.5 mag. 
In separating stellar objects from galaxies, we impose a rather strict quantitative criterion for 
the selection of stellar objects and separate stars and galaxies. To evaluate the objectiveness of 
the selection criterion, we simulate images with different sizes at different magnitudes and apply 
the selection criterion to estimate the fraction of stars classified as galaxies (incompleteness) and 
that of galaxies classified as stars (contamination). So the incompleteness and contamination 
are obtained as functions of image size and magnitude. Then we apply an incompleteness factor 
for each stellar object to estimate the true star count and apply a contamination factor for each 
galaxy to evaluate the number of contaminating galaxies. 



3.1. Selection Criterion 

We first simulate 500 PSFs with FWHM of 0."39 per 0.25 mag for J and 500 PSFs with 
FWHM of 0."29 for K' between 16 and 26 mag. Simulated objects are meant to be stars. Then 
we run the SExtractor to detect each artificial source and measure its magnitude, FWHM and 
stellarity. 

The details of the simulation procedure for each SExtractor run are as follows. We randomly 
place on the true object image 50 artificial objects with random counts within the magnitude 
range between m and m+0.25. We run the SExtractor using the same detection conditions as 
for the true objects and compare the resultant catalog with the true object catalog to create a 
list of new objects which newly appeared or brightened by more than 0.75 mag. We compare 
the coordinates of the 50 artificial objects placed prior to the SExtractor run with those of the 
resultant new objects and pick up the nearest neighbor pairs within 3 pixels (0."35) as recovered 
artificial objects. We repeat the same procedure ten times and move to the next magnitude bin. 
The number of the artificial objects, 50 is smaller than the total number of true objects, 350, and 
the confusion in the recovery process is due to the true objects. If 500 artificial objects were placed 
at once, confusion among them would dominate. We have avoided this situation by limiting the 
number to 50 for each SExtractor run, which is still much more efficient than the ideal case of one 
at a time. 

At J, the FWHMs of detected objects stay near 0."39 from 16 to 22 mag, grow gradually from 
22 to 23.5 and start to diverge at 24 mag. Below 25.5 mag, most of the sources are undetected. 
At K', the FWHMs of detected objects stay close to 0."29 from 16 to 21 mag, grow slowly from 
21 to 22.5 mag, and diverge at 23 mag. Below 24.5 mag, most of the sources are undetected. At 
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J, the stellarity stays close to unity from 16 to 23 mag, and then gradually drops below 0.8 at 
fainter magnitudes. Below 24.5 mag, only a few sources have the stellarity above 0.8. At K' , 
the stellarity stays close to unity from 16 to 22 mag, and gradually drops below 0.8 at fainter 
magnitudes. Below 24 mag, only one object has the stellarity above 0.8. From these simulations, 
we have found that J is more sensitive than K' for the sources with J — K' <1.0. In selecting 
stellar objects, we use the J band image, and the K' band image is used to measure the color. 

We further simulate 500 PSFs per 0.25 mag with FWHMs of 1.26xseeing, 1.26^xseeing, 
and 1.26^ X seeing between 16 and 26 mag for J. These sources are meant to simulate non-stellar 
objects. We again run SExtractor to detect each artificial source and measure its magnitude, 
FWHM and stellarity Simulated stellar objects are well separated from simulated non-stellar 
objects in FWHM down to 24 mag. Below J=24, stellar and non-stellar objects tend to mix 
toward fainter magnitudes. The stellarity stays above 0.8 at magnitudes brighter than J=17, 
18, and 21 respectively for PSFs with FWHMs of 1.26^ x, 1.26^ x and 1.26xseeing, and goes 
down near zero and comes back up to 0.7 below 24 mag. At brighter magnitudes, the FWHM of 
each artificial source is a good discriminator between a stellar and non-stellar object. At fainter 
magnitudes, the stellarity appears to be a superior discriminator because that of a non-stellar 
object rarely becomes above 0.7. 

From the simulations so far, we find following conditions for FWHM and stellarity to separate 
stellar objects from non-stellar objects for J. Stellar objects satisfy 

FWHM(") < 0.47 + loO-4(''-26) ^ (15) 

and 

Stellarity > 0.8. (16) 
The performance of the selection conditions is graphically shown in Figure 1. 

3.2. Selected Stellar Objects 

From the selection criterion given in the previous section, 14 objects are selected (Table 1). J 
and K' images of each object is shown in Figure 2. The completeness correction and the number 
of contaminants are estimated and given in Table 2. From Table 2, the star count is more than 
94% complete above J = 22.5 and 57% complete above J = 24. The completeness is only 28% 
between J = 24 and 24.5. After the completeness correction, the estimated total count is 24.3. 
The total number of contaminants is only 0.03 and is negligible. This indicates how strict the 
selection criterion is. The color-magnitude diagram of the stellar objects is given in Figure 3. 
Above J = 24, the J — K' colors of eleven objects are between 0.40 and 0.72. However below J = 
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24, one object is very blue (0.0) and two others are very red (> 1.0). 

4. Comparison of observations and models 

The completeness of the observed star count is 57% at J = 24. Model predictions and the 
completeness corrected counts are compared in Tables 3 and 4. The J — K' color of observed 
eleven objects brighter than J = 24 is confined in the range from 0.4 to 0.7. If the completeness 
correction is applied, the effective star count brighter than J = 24 is 13.5. 

The disk model A predicts only 2.6 stars brighter than J = 24, whose J — K color ranges 
from 0.78 to 0.96. We note that the color-magnitude relation is of the thin disk. 

The halo model B predicts 7.8 halo stars. Combined with the disk model A, the models by 
Gould ct al. based on HST observations predict 10.4 stars in the magnitude range between J = 18 
and 24, which agrees with the completeness corrected count at — 0.84(j. 

The disk model CI (thick disk scale height 910 pc) predicts 5.5 stars with J — K between 
0.87 and 0.96, while the disk model C2 (thick disk scale height 1.5 kpc) predicts 8.2 stars with 
J — K between 0.85 and 0.96. For both models, the thin disk contribution is small, 1.9 stars. 

The halo model D predicts 11.8 halo stars within the color range of 0.59 and 0.65. If some 
halo stars have colors of M subdwarfs instead of those of extreme M subdwarfs, the observed stars 
brighter than J = 24 may all be halo stars. In combining the disk model with the halo model D, 
the model CI is preferred because of the smaller predicted number count. The total predicted star 
count of Cl+D is 17.4 which is slightly greater than the completeness corrected count of 13.5 by 
1.0 a. The observed star count is between the prediction of models based on HST data by Gould 
et al. and that of models based on nearby stars LFs. Our limited statistics do not have strong 
constraint on the discrepancy between the results of two types of models. 

Apart from stars, the disk models predict a small number of substellar objects brighter than 
J = 24. The disk model A and C estimate the number of T dwarfs to be 1.4 and 0.39 respectively. 
Combining the local number density of T dwarfs with the Galactic structure of disk model A of 
Gould et al. apparently results in overestimating the number of T dwarfs. The disk model C is 
preferred for the combination with the local number density. The disk model A and C predict the 
number of L dwarfs to be 0.65 and 0.1. Since no object as red as L dwarfs brighter than J=24 was 
detected, the Galactic structure of the disk model C is again preferred. 



5. Individual Objects 
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5.1. Blue M subdwarfs, #7 and LHS 1826 

The object #7 is bluer {J — K' = 0.4) than any other object with J < 24. It is also bluer than 
any of the nearby extreme subdwarfs whose colors are given in Leggett et al. (1998)[ According to 



the calculations by Saumon et al. (1994)| of atmospheres of zero-metal low-mass objects, the J — K 



color of a zero-metal object near the stellar/substellar boundary will be blue and negative due to 
H2 collision induced absorption (CIA). The color of #7 may be affected by H2 CIA significantly. 
H2 CIA most likely affects the coolest and most-metal-poor object. In the solar neighborhood, the 
coolest and most-metal-poor object is LHS 1826 ( Gizis and Reid (1997)| ) whose infrared color was 



previously unknown. We obtained JHK photometry of LHS 1826 at the Okayama Astrophysical 



Observatory on 1999, November 30, using the facility infrared imager, OASIS (Okumura et al. 
(2000)1 ). Resultant magnitudes are J = 15.70 ±0.04, H = 15.43 ± 0.05, and K = 15.32 ± 0.07. 
Therefore J — K = 0.36 =b 0.08 and it is as blue as ^^7. It is important to measure the parallax of 
LHS 1826 and obtain its absolute magnitude. It may be the least luminous extreme M subdwarf 
known so far. ^7 and LHS 1826 indicate the presence of a faint blue halo population previously 
unnoticed. 



5.2. Red objects, #11 and #12 

Errors in J — K' colors of the objects with J > 24 are large. However, it is very likely that at 
least either #11 or #12 has J — K' > 1.0. Only L dwarfs are the stellar/substellar objects as red 
as these, but the expected number of L dwarfs between J = 24 and 25 predicted by the model CI 
is 0.06. 

Now we consider the possibility that at least one of the objects is extragalactic, namely an 
AGN. In the local universe, the fraction of AGNs to the entire galaxy population is estimated to 
be 1.3% based on the data of the CfA redshift survey ( puchra and Burg (1992)]) . There are some 
350 galaxies in SDF. If the fraction of AGNs does not evolve, there may be four AGNs in the field. 
However, we know that galaxies evolve and so do quasars (e.g. ISchmidt et al. (1995)| ). So the 



validity of the number estimate based on the local universe information is questionable, 

PG quasars are red and most of them have J — K between 1 and 2 ( [Neugebauer et al 



(1987)1 ). In terms of colors, low-redshift quasars are candidates. However they are extremely 



bright (J <12). If we assume that the object has a power-law continuum whose spectral index 
a = —0.653 is consistent with J — K' = 1.5, the combinations of absolute B magnitude, Mb and 
z can be obtained for the apparent magnitude of J = 24.2. Here the cosmological parameters we 
assume are Hq = 70 kms~^ Mpc~^, = 0.3, and Vl\ = 0.7. Allowed combinations of {Mb,z) are 
(-22,1.2), (-24,2.9), and (-26,7.0). So if the object is a quasar with Mb = -26, it must be located at 
very high redshift. If such quasars are as numerous as to be detected in a small field, the number 
density may be high enough to affect reionization of the universe. 
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6. The Galaxy Seen by NGST 

Since we have obtained the star count models which are capable of predicting the appearance 
of the Galaxy at faint magnitudes, it is interesting to use them to envision the Galaxy seen by 
NGST. In Table 5, surface number densities (arcmin~^) of objects between K = 27 and 30 are 
given for the north Galactic pole. 

In the galactic disk, T dwarfs are by far the most dominant population. However the expected 
number count is small (< 0.01) due to the small vertical scale height. There may be fainter and 
bluer M subdwarfs than predicted, because the LFs in the models extend only to My = 14, while 
those of the disk extend to My = 18. Old halo white dwarfs are the most numerous of all stellar 
populations. T dwarfs are blue in the near infrared due to molecular absorption bands. Cool M 
subdwarfs and halo white dwarfs are blue due to H2 CIA. Stellar astronomy will have a blue faint 
end in the near infrared. 

7. Conclusions 

Our primary conclusions follow. 

(1) We constructed infrared star count models of the Galaxy applicable at faint magnitudes 
based on the models and luminosity functions obtained at V. There are two types of models 
derived in our paper. First we derived infrared versions of the disk and halo models by Gould et al. 
who obtained the local luminosity functions and functional forms of the models based on HST star 
counts. Second we constructed a double exponential disk model based on the luminosity function 
of the nearby stars by Reid and Gizis and a halo model based on the nearby subdwarf luminosity 
function of Dahn et al.. In addition to dwarfs and subdwarfs included in original models, we also 
took into account L dwarfs, T dwarfs, and white dwarfs of both disk and halo. 

(2) As a test case of the models, we analyzed the infrared imaging data at J and K' obtained 
during the Subaru Deep Field (SDF) survey to study stellar objects. Out of some 350 objects, 
14 stellar objects were selected in a 2' x 2' field based on a morphological criterion applied to 
the J band image. Both completeness and contamination associated with the selection criterion 
were evaluated by simulations. The J band image is 57% complete at J = 24 and the number of 
contaminants was estimated to be negligible. The prediction of the HST based models agrees with 
the observed count at — O.Scr and that of the nearby star LF based models also agrees with the 
observations at +1.0(7. The observed count is between the predictions of the two types of models 
that have contradictory local luminosity functions. With our limited statistics, the observing data 
do not favor a particular type of models. 

(3) Infrared star count models we obtained predict that the Next Generation Space Telescope 
will primarily see T dwarfs, M subdwarfs, and old halo white dwarfs at faint magnitudes. 
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Table 1. Stellar Objects 



# 


(a,5)(J2000) 


J 


a(J) 


FWHM (J) 


Stellarity 


J-K' 


a{J-~K') 


1 


132229.6+272841.3 


21.999 


0.016 


0.420 


0.990 


0.518 


0.04 


2 


132240.0+272932.1 


18.217 


0.001 


0.435 


0.980 


0.696 


0.0014 


3 


132250.6+272900.9 


20.209 


0.003 


0.418 


1.000 


0.714 


0.006 


4 


132304.3+272904.3 


20.269 


0.004 


0.416 


0.980 


0.523 


0.008 


5 


132309.5+272949.8 


21.389 


0.009 


0.435 


0.980 


0.678 


0.014 


6 


132314.1+272941.9 


20.869 


0.006 


0.421 


0.980 


0.635 


0.0116 


7 


132315.8+272921.2 


23.545 


0.048 


0.416 


0.960 


0.403 


0.10 


8 


132344.9+272841.3 


22.191 


0.021 


0.461 


0.970 


0.669 


0.039 


9 


132357.7+272951.8 


21.963 


0.014 


0.419 


0.980 


0.473 


0.030 


10 


132525.4+272922.4 


23.886 


0.057 


0.386 


0.970 


0.661 


0.115 


11 


132558.9+272854.2 


24.215 


0.121 


0.423 


0.890 


1.288 


0.17 


12 


132607.4+272907.0 


24.433 


0.115 


0.474 


0.870 


1.648 


0.16 


13 


132617.2+272945.4 


23.804 


0.066 


0.501 


0.950 


0.721 


0.141 


14 


132622.1+273010.3 


24.176 


0.109 


0.458 


0.810 


0.002 


0.255 
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Table 2. Completeness correction and contamination 



J 


Detection 


Corrected Detection 


Contamination 


17.2500 


0.0000 


0.0000 


0.0000 


17.7500 


0.0000 


0.0000 


0.0000 


18.2500 


1.0000 


1.0142 


0.0000 


18.7500 


0.0000 


0.0000 


0.0000 


19.2500 


0.0000 


0.0000 


0.0000 


19.7500 


0.0000 


0.0000 


0.0000 


20.2500 


2.0000 


2.0555 


0.0000 


20.7500 


1.0000 


1.0277 


0.0000 


21.2500 


1.0000 


1.0320 


0.0000 


21.7500 


2.0000 


2.0725 


0.0000 


22.2500 


1.0000 


1.0638 


0.0000 


22.7500 


0.0000 


0.0000 


0.0000 


23.2500 


0.0000 


0.0000 


0.0020 


23.7500 


3.0000 


5.2724 


0.0180 


24.2500 


3.0000 


10.7527 


0.0100 


2i.7r)()() 


().()()()() 


O.OOIK) 


I). ()()()() 
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J 


disk model A 


halo model B 


K+B 


Corrected count 


18-20 


1.53 


0.28 


1.82 


l.Olil.OO 


20-22 


0.85 


1.86 


2.72 


6.19±2.49 


22-24 


0.23 


5.68 


5.91 


6.34±2.51 


Total 


2.61 


7.82 


10.45(-0.84c7) 


13.54±3.68 



Table 3: Models by Gould et al. 
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J 


disk model CI 


halo model D 


Cl+D 


Corrected count 


18-20 


2.13 


0.40 


2.74 


l.Olil.OO 


20-22 


2.32 


3.00 


5.31 


6.19±2.49 


22-24 


1.08 


8.43 


9.50 


6.34±2.51 


Total 


5.53 


11.83 


17.36(1.04a) 


13.54±3.68 



Table 4: Models based on nearby star LFs 



-17- 



Stellar population disk model A 


halo model B 


disk model C halo model D 


T dwarfs 0.002 




0.005 


M subdwarfs 


0.50 


0.12 


Halo white dwarfs 


1.53 


1.23 



Table 5: Predicted surface number density (arcmin ^) 
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Fig. 1. — Selection of stellar objects. In the upper panel, small dots indicate simulated objects 
for FWHMs 0."39, 1.26 x 0."39, (1.26)2 x 0."39, and (1.26)3 x 0."39. 500 objects are simulated for 
each 0.25 mag interval. Dots generated for different FWHMs tend to mix at fainter magnitudes. 
Observed objects with different J magnitudes are indicated by circles with various sizes. The 
selection criterion for FWHM is shown by the solid line below which the sources are regarded 
as stellar objects. In the lower panel, the stellarity condition is imposed on both observed and 
simulated objects and only those that satisfy the condition are plotted. When both the FWHM 
and stellarity conditions are imposed, the mixing of dots generated for different FWHMs is fully 
suppressed. 

Fig. 2. — J and K' images of stellar objects (The J image is on top of the K' image in each pair). 
Pairs of ^^1 through ^5 are presented in the bottom raw from left to right, those of ^7^6 through 
#10 in the mid raw and those of #11 through #14 in the top raw. 

Fig. 3. — J YS J — K' diagram. All the eleven objects with J <24 have J — K' colors in the range 
of 0.4 and 0.7. All the objects with J > 24 have exotic colors. Error bars are of 1 a. 
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